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Abstract 


This report summarizes research cmidiic led during I he past three and a half 
years aimed al developing and testing a I urhulence/l ransition model applicable 
to high-speed turbulent flows The lirsl two years of the project focused on 
fully turbulent Hows, while emphasis shifted to boundary layer development in 
the transition region during the final year and a half. 'This report includes a 
brief summary of research accomplished during the iirst three years and cites 
publications that describe research results in greater detail. The main body 
of this report summarizes research conducted during the final six months of 
the period of performance. The primary results of the last six months of the 
project are elimination of the H'-w model's sensitivity to the freeslreum value of 
lj and development of a method for triggering transition at a specified location, 
independent of the freest ream turbulence level. 
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Chapter 1 


Introduction 


The purpose of t his research project Tor the past three and a hah years has been 
to develop analytical and computational tools suitable Tor predicting propert.es 
of hypersonic flows, including both turbulent and transitional flow regm.es. Our 
research efforts have generated a total or eight publications, including this final 
report. The previous publications are as follows. 


I Wilcox, I). “Hypersonic Turbulence Modeling Without the Kpstlon 
Equation,” T h National Aero-Space Plane Technology Symposium, Ses- 
sion VII, Paper No 27 (October 108!)) 

2. Wilcox, I). 0., “A Half Century Historical Heview of the k-u> Model," 
AIAA Paper 91-0015 (January IIMH) 


3. 


Wilcox I) 0 “Theoretical Study of Turbulent Mixing Ihtween Hyper- 
sonic Streams’” DCW Industries It. port No DOW- IT-35 01 (May 15)91). 


i Wilcox, I). O., “Progress in Hypersonic Turbulence Modeling, AIAA I a- 
per 91-1785 (June 1991) 

5 Wilcox 1). 0. “Progress in Developing a I ransition Model Tor Iligli-Spn.d 

Flows,” DOW Industries Report No. DOW IT-35-02 (October 1992). 


G. Wilcox, D C., “'The Remarkable Ability „r Turbulence Model Equations to 
Describe Transit ion,” Fifth Symposium on A' win rical and I'hysuul Aspids 
of Atwdynami, Flows, Oalilbrt.ia State University, bong Head., California 
(13-15 January 1992) (accepted for publication in the AIAA Journal). 

7 Wilcox DC., “Dilatation Dissipation Corrections Tor Advanced Ihrbu- 

Icnce Models," AIAA Journal, Vnl 30, No. II, PP 2639 2040 (November 
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Publications I, 3, 4 ami 7 deal primarily with fully-lurbulenl, high-speed 
Hows, and represent research results Tor the first two years of the project. One 
of the primary accomplishments was discovery that separating dissipation into 
soleuoidal and dilatation contributions improves predictive accuracy for the 
compressible mixing layer at the expense of a loss in predictive accuracy for 
boundary layers. We identified the cause of this problem and postulated a 
dilatation-dissipation correction that is accurate for the mixing layer and the. 
boundary layer. A second accomplishment was identificat ion that the k- v model 
faithfully reproduces the compressible law of the wall for high-speed flows while 
the k-i model does not Perhaps the most important accomplishment was the 
demonstration, in Publications 1 and 7 that dilatation dissipation does little to 
improve two-e<juatiou-liiodel predictions for shock-separated Hows. Rather, the 
research verilies that using a second order closure model greatly improves pre- 
dicted size of the separation region, although oilier How properties are not accu- 
rately predicted (e g., velocity profiles downstream of the rcaUachment point). 

Publications 2, \ and (i introduce low-Ueyuolds-nmnber corrections for the 
k-ut model that permit the model to accurately describe near-wall properties 
of high-speed boundary layers. The modifications presented in the first two 
publications ultimately have been superseded by those postulated in Publication 
G. The latter Publication represents much of the btisic research conducted during 
the third year of the project. 

Publications f> and G focus on application of the low-Reynolds-uumber k-u> 
model to a collection of transitional Hows. The publications show that the 
model accurately predicts properties in the transitional region provided tran- 
sition Reynolds number is relatively small iiy contrast, predicted transition 
width is generally smaller than measured for transition at high Reynolds num- 
bers. Additionally, computed results show a stronger-t han-desircd sensitivity 
to the Ireestream value of w, and the need for a simple method of triggering 
transition at a specified location. 

During the final six months of the project, we have eliminated two of the 
three outstanding problems with tin* model as it applies to transitional bound- 
ary layers. Specifically, we have devised a modilical ion that eliminates solution 
sensitivity to the freesiream value of uj , and we have devised a straightfor- 
ward method for triggering transition at a specified location. However, limited 
funding has precluded solving the problem of the too-ahrupt transition at high 
HrynohiH numbers 

As summarized in Publication f>, the overall objective of our approach to 
the transit ion problem is to use tin Wilcox 1 k -u> turbulence model as the foun- 
dation for studying and modeling the transitional flow region. Consistent with 
the needs of NASA, the transit ion point is assumed to be known a priori. Com- 
pulations can thus be initiated at the known transition location and continued 
downstream through the transitional How region and well into the fully turbu- 
lent region To develop the model, we have followed a sequence of interrelated 
steps. 





1. Ill Publications 2, A ami (i, wo drew from (hr extensive work dour by 
Wilcox et al. 2 “ 7 to help formulate a low-Reynolds-number version of the 
k-w model. 

2. In Publications A ami f>, we simulated several of the Hows that have been 
done with Direct Numerical Simulation (DNS) methods. I he results oT 
Mansour, Kim and Moin 8 were computed for two different Reynolds num- 
bers. Comparing model predictions with the DNS results helped greatly 
in developing the model. 

3. In Publication 5, we tested the model against It) two-dimensional, fully 
turbulent boundary layers, including both incompressible and compress- 
ible cases. We have repeated 5 of these cases in this report. 

A. In Publication 5, we tested the model against all or the two-dimensional 
cases presented by Singer, Dinavahi, and Iyer. 0 We have repeated ail of 
these cases in this report. 

5. We have analyzed the compressible log layer to explain why adding a cross- 
diffusion term to the u> equation can have an adverse effect on predicted 
skin friction. 

As noted above, addressing Items I and ‘2 lias produced a low-Heyuolds- 
iiuinber version of the Jt-ui model. Publication 0 describes the model and 
presents applications to fully turbulent channel and pipe flow, and for a tran- 
sitional, incompressible flat-plate boundary layer In this report we propose a 
modified version of the model that eliminates freestream boundary-condition 
sensitivity. This report applies (he mode) to turbulent boundary-layer test 
cases, two free shear flow cases and to more than 20 transitional Hows for which 
experimental data are available 

Chapter 2 summarizes the revised low-Reynolds-number version of the k-w 
model. Chapters 3, A and f> present results of our applications. Chapter fi 
summarizes results of the research. The effect of a model revision on the com- 
pressible law of the wall is discussed in Appendix A. 
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Chapter 2 


Revised Low Reynolds 
Number k-w Model 


The revised low Reynolds- number version of I lie k-u/ model is wrillen in terms 
of l'avre l(J mass averaged variables. The equat ions of motion are as follows. 
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In Initiations (2.1) (2 5), f and x x denote lime ami position vector; p and 
p are density and pressure; a, is mass averaged velocity vector; E and U are 
total energy and total enthalpy; r tJ and </, arc the total stress tensor and the 
total heat flux vector; is the Meynolds stress tensor; k and u are turhulence 




kinetic energy and specific dissipation rale; and, a, /?, /i*, *t, rr* , are closure 
coefficients. The following constitutive relations are needed to close the system. 
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In the constitutive relations, we have introduced the molecular stress tensor, 
, mean strain-rate tensor, ,S’ f; , mnleriilar viscosity, /i, eddy viscosity, p?-, 
enthalpy, /i, laminar Prandtl number, /V/ , turbulent Prandtl number, /Vr, 
internal energy, e, and an additional closure coefficient, u* . Finally, to complete 
closure of the system, we specify the values of the closure coefficients »is follows. 
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The (plant ity llvj' is the turbulence Reynolds number defined by 
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J lie primary difference between tins model and llie model introduced by 
Wilcox 11 is the introduction of the term in liquation (2.5) proportional to cr d . 
1 his term is similar to Menter’s 1 " cross diffusion term. The effect of this term is 
to replace the entrainment velocity, iq in the equation by (v ~ Ok/dy). 

Since k decreases approaching the shear layer edge (assuming a d > 0), the 
net ellect is to make the effective entrainment velocity positive (or at least less 
negative). As a result, u> diffuses from the turbulent region into the nonturbulenl 
region, which is the opposite of what happens with the k-u> model. Thus, the 
freest real u value of u has no effect on tin* solution 

As shown by Wilcox 13 , it is important to suppress the cross -diffusion term 
close to solid boundaries lor wall bounded Hows. The prescription in Equa- 
tion (2.15) causes Od to vanish near a solid boundary since k increases and u 
decreases in the viscous sublayer. This* mociilication to the u> equation eliminates 
the model s sensitivity to the Irei sl ream value oi w. The values of several clo- 
sure coefficients must be modified to achieve optimum results fur both boundary 
layers and free shear layers. Specifically, a and a* must assume values larger 
than used in the Standard k-*j model This, in turn, requires changing a and 
to maintain a satisfactory law of the wall in the turbulent boundary layer. 
The value quoted for a guarantees that the K dr man constant is 0.41 while the 
value ol lt w yields an additive constant in the law of the wall of 5.0. In summary, 
Jable 2.1 shows how the revised models closure coefficients differ from those of 
the Wilcox 11 model. 


Table 2.1: Closure Coefficient Differences 


( Coefficient 

New Model 

Wilcox 11 Modd 

a 

3/5 

1/2 

< 7 * 

1 

1/2 

Od 

Equal ion (2.15) 

0 

Or 

1/2 

r./a 

IL 

N/5 

27/10 
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Chapter 3 

Turbulent Boundary-Layer 
Applications 


A key objective of research during I fit' filial year and a half of this project has 
been to describe boundary- layer development t hrough transition from laminar 
to turbulent flow. Of course, it is important that we provide an accurate de- 
scription in the turbulent region immediately following transition, i.e., we insist 
that our model approach the proper limiting state of the boundary layer. Con- 
sequently, since the k-w model without viscous corrections generally is very 
accurate for turbulent boundary layers, a round of tests is in order to make sure 
model predictions are not adversely affected by the viscous corrections. VVe 
have performed five boundary layer compulations including effects of adverse 
and favorable pressure gradient, and lor Mat h numbers as high as 10. Table 3.1 
lists the five cases. 

I able 3.1: Turbulent Boundary Layer Test Cases 


Description 

V,, 

Data Source 

Incompressible Flat Plate 

None 

YVieghardt-Ti!lmaii M 

Samuel Joubert 

A d verse 

Samuel Joubert 15 

Much 2.244 Flat Plato 

None 

Shulls'® 

Mach '1.544 F lat Plate 

None 

Coirs' 6 

Macl. 10.31 Flat Plalo 

None 

Walson 16 


Incompressible Flat Plate Figure 3.1 compares computed and measured 
skin friction, cj , and sublayer scaled velocity, u f As shown, the only significant 
difference between predictions with and without viscous corrections is in the 
transition point I he unmodified model undergoes t ransition much earlier than 
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the model with viscous corrections. Skin (rid ion at the final station is 2.51 10~ 3 
witli viscous modifications compared to 2 1 ( J 10~ 3 without. 

Samuel- Joubort. This incompressible adverse pressure gradient boundary 
layer was supposed to be a simpler application in Stanford Olympics II. 17 On 
the contrary, it proved to be very difficult for all models and has become a key 
lest case for bow well a turbulence model predicts effects of adverse pressure 
gradient, f igure 3.2 shows that both high and low-Heynolds-nuinber versions 
of the model nearly duplicate measured skin friction and velocity profile at the 
last station. At the final station, cj = 1.35 IIT 3 with viscous corrections and 
1.25 I(T 3 without. 

Compressible Flat Flutes As shown iu Figures 3.3, 3.4 and 3.5, the 
viscous corrections have virtually no died on skin friction aud velocity profiles 
in the fully turbulent region for these three applications. For all three cases, 
skin friction differs by less than one percent at the final station. 

In summary, for all 5 cases considered, dilferenres between computed flow 
properties with and without viscous corrections are less Ilian 8%, ami are gen- 
erally less than 2%. Additionally, although not shown here, we have computed 
(under separate funding) I f more incompressible boundary layers, 3 with fa- 
vorable pressure gradient and 11 with adverse pressure gradient. Along with 
the (lal-plate boundary layer and the Samuel Jouberl case, these are the 1G 
boundary layers analyzed by Wilcox 1 *. As with the cases discussed above, the 
low-Hcynolds-number model is ;is accurate as the basic model without vis- 
cous modifications. Hence, the viscous corrections leave the best features ol 
the F-w model intact, i.e., the model still accurately predicts elfects of pressure 
gradient and compressibility up to Mach 5. 

While the cross-diffusion term in the u> equation greatly reduces the model’s 
sensitivity to the freest ream value of cj, it does not entirely remove it. A larger 
value of ad is needed to complete ly remove the sensitivity. However, as shown 
in Appendix A, using a larger value for a,t yields an inaccurate prediction for 
the compressible law of the wall. 
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Chapter 4 

Free Shear Flow 
Applications 


While k - w model boundary-layer computat ions display a mild sensitivity to the 
freeslream value of w, the sensitivity is much stronger for free shear flows. As 
shown by Wilcox 1 **, the spreading rale lor I he incompressible mixing layer varies 
between 0.103 and 0.141; the measured spreading rate is 0.115 Similarly, for 
the far wake, computed spreading rate varies from 0 301 to 0.500, as compared 
to a measured value of 0.305. 

i he new model almost completely eliminates this problem. 1'igure 4.1 com- 
pares computed and measured mixing-layer velocity profiles lor two values of 
the Irecstream uj that differ by several orders of magnitude. The computed 
spreading rate varies between 0.105 and 0 107. 

Similarly, figure 4.2 compares computed and measured velocity profiles for 
the far wake. For a wide range ol freeslream values of u>, the spreading rate 
varies bet ween 0.285 and 0.318. 

We have attempted compulations for the plane and round jets, but have 
encountered numerical difficulties in obtaining converged solutions. This illus- 
trates a key problem attending use of the cross-diffusion term. Jn general, the 
term greatly increases the stiffness of the model equations by altering the effec- 
tive time scale 1 ) near a turbulent/ uoul iirbulent interface. As a consequence, 
the convergence rate of the numerical algorithm tends to be reduced. 











Chapter 5 


Transitional 

Boundary-Layer 

Applications 

5.1 Triggering Transition 

I he question of how to trigger transition for a given application raises an impor- 
tant. point about using turbulence model equations Tor transitional flows. YVe 
must keep in mind that transition is a complicated phenomenon 'IVansition 
* s by ft disturbance in a boundary layer only if the frequency of the 

disturbance falls in a specific hand Reynolds averaging has masked all spectral 
effects, and all the model ran represent with k and lj is the intensity of the 
disturbance and an average frequency. Hence, it is possible for the turbulence 
model to predict transition when it shouldn't occur. If we choose to trigger 
transition via the presence of a freest ream disturbance, using turbulence-model 
equations js sensible only if the triggering disturbance is broad band, i.e., con- 
tains all frequencies. 

In devising the model, the closure coefficients «• and a 0 [Equation (2. 16)] 
have been chosen to match the minimum critical Reynolds number for the Bla- 
sius boundary layer Consequently, the model can be expected to yield sensi- 
ble transition predictions for constant pressure, incompressible boundary layers. 
Unfortunately, the model equations must be recalibrated (i.e., o* and o fJ must 
be adjusted) to accommodate each new complicating effect [sec YVilcox el al, 2 “ 7 J. 

In this project, the goal has hern lo assume the transition point is given and 
to use the model equations t hrough the transitional region and into the turbulent 
part of the flow. To date, the freest ream turbulence intensity, proportional to 
the freest ream value of A’, has been adjusted to match the measured transition 


o* r 


18 



point. This is satisfactory when the* transition point occurs at a large Reynolds 
number, which requires k^ to be small relative to However, for high-speed 
Hows in which transition occurs at a relatively small Reynolds number, we have 
found that unreasonably large values of k iXj are needed to trigger transition, so 
large as to afTect the total energy in the freestreain in a physically unrealistic 
manner. Thus, a new method for triggering transition is needed. 

To devise an alternative method, we can take advantage of a unique feature 
of the k-LJ model. Specifically, by using a finite value for u) at the surface, we 
can simulate surface roughness with the model. For fully turbulent boundary 
layers, Wilcox 1 shows that 

w = — Sh al y = U (5.1) 

v 

where u r is friction velocity and Sn is a dimensionless function of the surface 
roughness height, kn t defined by (with kj^ =. Ujkn/v): 

[ (BU /kj t )\ *£<25 

.V* - l (5 2) 

[ 1 Oll/ir/ii *£ > 25 

Since increasing the surface roughness height corresponds to decreasing the sur- 
face value of u (and thus the dissipation in the k equation), the model predicts 
that roughness will have a destabilizing eU’cct. This is consistent with measure- 
ments, and patches of surface roughness are often used to trigger transition in 
experiments. Thus, a possible way to trigger transition with the model equa- 
tions is to numerically simulate a roughness strip via Equations (5.1) and (5.2). 

We have run more than 2(1 two-dimensional transitional boundary layer cases 
to test this idea; results of the computations arc given in the next section. We 
have been able to trigger transition at the desired location for all of the cases 
considered using a roughness strip with k^ and the slreamwise extent of the 
strip, As, given by (lie following correlations. 



The quantities b t and ltt Si are the boundary-layer thickness and transition 
Reynolds number based on arc length. 


5.2 Applications 

In order to test the new transit ion- triggering method, we have computed all of 
the two-dimensional transitional boundary layer cases considered by Singer, el 


hI. Our applications also include a MjicIi 20 rf'-< i nlry raso. fable 5.1 summarizes 
the rases we have done. 


fable 5.1: I ransilional Boundary Layer lest Cases 


f Flow 

Description 

Vp 

Data Source 

1 

Incompressible Flat Plate 

None 

Scbiibaiier K lebanoff 1 9 

2 

Favorable Vp 

Fa vorable 

Wair-W»*rli* 20 

3 

Supersonic Cone Flow 

None 

Fisher- Dougherty 21 

5 

Supersonic Cone How 

None 

< -hen, rt a I. 22 

4 

Freestream Turbulence 

None 

Scbiibaiier Sk rams tad 2 * 1 

4 

Freestream 'hirbiilence 

None 

Blair 2/1 

5 

Prolate Spheroid 

favorable 

Meier- Kreplin Ming 25 

6 

Concave Surface 

None 

Swearingen Black welder 20 

G 

Mach 20, 5 f> Cone 

None 

Howard 27 


In all cases, computation begins at I lie* plate leading edge, and the turbu- 
lence kinetic energy is initially set to an extremely small value, viz., I0" ,fi F/£, t 
throughout the boundary layer. This value is too small to trigger transition 
naturally. I he initial u> profile is given by tin* exact laminar-flow solution to the 
model equations. 1 1 

In addition to testing the transition triggering method described above, we 
have tested the effect of the cross-diffusion term in Equation (2.5). A disap- 
pointing feature of our previous transition predictions was a mild sensitivity of 
transition location to the initial w profile. The sensitivity was caused by the 
small, but noticeable, effect of t he freestream value of u ; on the laminar-flow so- 
lution for u;. I fie primary reason for introducing the cross-diffusion term was to 
eliminate this sensitivity for transition predictions. As shown in the preceding 
chapters, the cross-diffusion term eliminates most of the sensitivity for turbu- 
lent shear flows. Numerical experimentation For several of the transitional cases 
listed in Table 5.1 verifies that this is true for transitional flows also. That is, 
our computations verify that transition predictions are completely insensitive to 
the initial u> profile, even when transition is triggered by freestream turbulence. 

Incompressible Flat Plate The first case is Flow I from Singer, et al. 9 
I Ills is an incompressible flat plate boundary layer that undergoes transition 
at a plate length between 1.0 in. and 18 m According to Fq nations (5.5) and 
(5.4), this How requires n roughness strip with kn fh x — 5.0 and A s/S t = 25. 
Numerical experimentation shows that a shorter transition strip, i.e., A.s/A* = 7 
is sufficient to trigger transition at the desired location for this case. 

figure 5.1 coinpan's computed and measured skin friction throughout the 
transition region for the model with and without the cross-diffusion term. In 
both cases, computed and measured cj differ by less than 10% of the peak skin 
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I igure 5.1: lVaitbit ioual incompressible Hut-plate boundary layer; New k-u> 

model; k-u) model witboul cross diifusion. 

Iri< lion. Note that with cross -dill usion included, the distance between minimum 
and maximum skin friction is a bit less than that without cross diffusion. 

l*nvornl>lo Pressure Gradient lloundary Layers The next applications 
art' for incompressible boundary layers in a lavorable pressure gradient The 
boundary layers considered correspond to Flow 2/Gases 1 and 2 of Singer, ct al. 
In addition to having adverse pressure gradient, t lie surface is cooled. Figure 5.2 
compares computed and measured Stanton number for the two cases. (last* 2 
has a stronger favorable pressure gradient than Gn.se 1 The dimensions of the 
roughness strip required to match the measured transition point are (knfb tt 
*\fi/b t ) — (10, 10) and (1, 0), respectively, for ('uses I and 2. Filiations (5.3) 
and (5/1) indicate (kn/b lt Aa/^ ( ) = (8.5, 8.7) and (6.0, 12.5), respectively. 
Again, we see that the transition width is reduced a bit when the cross-diffusion 
term is included 

Supersonic Cone Plow, for this application we consider Flow 3 of the 
Singer, el al. study. Gases I I focus on How past a 5° half angle cone with 
Mach numbers ranging from l.ili to 1.86, corresponding to measurements of 
Fisher and Dougherty.- 1 < Ws 5 7 are for Mach 3,30 llow past a 5° half-angle 
cone , with an adiabatic surface (see (’hen, et ul. J “). Table 5.2 summarizes the 
transition Reynolds number, /(e X|l and Reynolds number bused on transition 
width, ltc±, t lor all 7 c«ises, lor the new and old low Reynolds-uumber k-w 
models. 
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Table 5.2: Results for Fisher- Dougherty and Chen el al. 'lest Cases 


| Case 

Mach Number 

Ite 

il 

(IUSj 

't ) nr m 

(HtAx t ) v i4 

(Mr A* 


mm 

ltd 

7.13 

10 6 

0 45 

10 6 

0 89 

10“ 

069 

10“ 

ii 

1.30 

5.70 

10“ 

0.50 

10“ 

0.84 

10“ 

0.97 

10“ 

■1 

1 55 

7.00 

10“ 

0.72 

10“ 

1 .08 

10“ 

1 24 

10“ 

mm 

1.80 

7.49 • 

10“ 

0.78 

10“ 

1 .07 

10“ 

1 56 • 

10“ 

HI 

3 30 

6.92 

10° 

1 01 

10“ 

I 63 

10“ 

— 


u 

3.30 

9.51 • 

10“ 

1.19 

10“ 

1 79 

10“ 

— 


m 

3.30 

1.07 

I0 7 

1 30 

10“ 

1 .82 

10“ 




Finest ream Turbulence Effects. We turn now to effects of freestream 
turbulence, i.e, I low 4 of the Singer, el al. study. AH of the cases considered are 
for incompressible boundary layers, 'table 5.3 shows that computed transition 
widLh for Cases 1 5 is consistently smaller than measured, and noticeably smaller 
than the model without cross diffusion. 

Table 5.3: Results for Schubauer-Skramslad 'lest Cases 


Case 

V 

i Rex, 

{Ur 

f )n r II' 

(Rcaj 

, )»M 

(Kt&j 


1 

.012% 

2 85 

10 e 

0.56 

10 s 

0.88 

10 e 

1.00 

10® 

2 

.100%. 

2.75 

10“ 

0.52 

10“ 

0.86 

10“ 

1 .20 

10® 

3 

.200%. 

2.20 

10“ 

0 18 

10“ 

0.78 

10“ 

1.50 

10® 

4 

.260%. 

1.80 

10“ 

0 15 

10“ 

0.69 

10“ 

1.40 

10® 

5 

.310%, 

1.40 

10“ 

0.36 

10“ 

0.58 

10® 

1.25 

10® 


For Cases (i, 7 and 8, Figure 5.3 shows that predicted peak skin friction is 
about 10%-15% lower than measured. The results for these three cases are very 
similar to those obtained without cross diffusion. 

Incompressible Flow Past a Spheroid This case is Flow 5 from Singer, 
et al. The flow examines transition of a boundary layer over a prolate spheroid 
at zero angle of attack, with transition triggered by a roughness strip. Ex- 
perimental data have been provided by Meier, Krcpiin and Ming. 25 Table 5.4 
summarizes the four cases. 

Figures 5.4 through 5.7 compare computed and measured skin friction (based 
oil local bou udary layer edge velocity) ;is a function of arc length along the 
spheroid. As shown, computed and measured skin friction are closest when 
transition is induced by the roughness strip. Without the rouglmess strip, the 
predicted transition occurs more abruptly than measured. 


















iable 5.4; Meier Kroplin Ming Test Cases 


C a-se 

Description 

kh/ht 

A a/bt 

1 

Uckj = 20 m/sec, no roughness strip 


6.0 

2 

= 20 m/sec, with roughness strip 


8.0 

3 

= 30 m/sec, no roughness strip 


1.6 

4 

— 30 m/sec, w it h roughness st rip 


ISO 


Concave: Surface Boundary Layer In this application we consider in- 
compressible How over a concave surface. \Ve present results in Figure 5.8. As in 
all previous applications, including cross diffusion causes a noticeable reduction 
in transition width According (o Filiations (5.3) and (5.4), the transition strip 
should have kn/t t — 3.5 and A a/6 t = 21. Nuinericai experimentation shows 
that somewhat smaller values, k H /b t = 3 0 and A*/d, = 12, can be used. 

Mach 20, 5° Corn*. Figure 5.0 compares computed and measured Stan- 
ton number, Si t lor Macli 20 llow^ 7 past a 5° ball-angle cone. As shown, the 
computed Stanton number increases more abruptly than measured, although 
less abruptly with cross dillusion included. This is the only case for which 
transition width increases when cross diffusion is included. Note also, that the 
computed Stanton number downstream o! transition is closer to measurements 
with cross dillusion included. FijuaLions (5.3) and (5.4) indicate a transition 
strip with kn/d t = 3 0 and As/6, = 152 Numerical experimentation shows 
that transition can be triggered with k fi /6 t = 1.0 and A s/b t = 32. 



25 



















Chapter 6 

Summary and Conclusions 


lu tlie context o| turbulent boundary layers, ibis report shows that the new 
low- Heynolds- number version of the model retains most of the best features of 
the baseline, high Heynolds-miiuber version of the model. Most notably, the 
model is just as accurate us the baseline model for Hows with adverse pressure 
gradient, i his guarantees that the model approaches the correct post- transition 
asymptotic state. 


1 h<r cross-dill usion term is very effective in eliminating solution sensitivity 
to the freeslreum value of w for turbulent boundary layers, free shear flows and 
transitional boundary layers. However, (his has been accomplished at the cost 
of a noticeable increase in the stillness of the equations of motion. 

While the new transition- triggering method has proved to be very effective, 
and the cross-diffusion term eliminates transition sensitivity to the initial u/ 
profile, computed transition width is actually reduced somewhat with cross 
diffusion included liquation (5.3) is an upper bound for the roughness height 
required to trigger transition for all of the cases considered. In many cases, a 
smaller roughness height is sufficient, and the user should try different roughness 
heights, if possible, to determine the optimum height for a given application, 
liquation (5.4) is less certain The formula also expresses an upper bound that 
covers all ol the cases considered. I he actual values used permits the roughness 
strip to persist for at least three slreamwjsc linite-differcnce cells. 



Appendix A 


Compressible Law of the 
Wall 


While the cross-diffusion term is very effective in eliminating the k-u model’s 
sensitivity to the freostream value of w, it has a negative impact on predictions 
for compressible boundary layers. Figure A.! compares computed and measured 
velocity profiles for Mach 1.5 and Mach 10.3 flat-plate boundary layers. The* 
computations have been done with all closure coefficients chosen to insure that 
the law of the wall holds in the limit of incompressible flow, viz, 

ti + = + fl ( A. 1 ) 

with k = 0.41 and H — 5.0. As shown, when cross diffusion is omitted (cr rf = 0), 
the law of the wall (in terms of Van Driest 38 variables) is reproduced, and 
close agreement with measured velocities is obtained. By contrast, when <r<t 
is 0.5 (with (Td/cr — 1), significant differences between computed and measured 
velocities appear at both Mach numbers. The figures indicate smaller differences 
for <Td = 0.3 (with rx rf /<r = I / 2 ) , and this is the value used in all computations 
of Chapter s 3 through 5. 

To explain the adverse effect of cross diffusion on the compressible law of 
the wall, this appendix examines k-u> model predicted compressible log-layer 
structure. The results are particularly illuminating and clearly demonstrate 
why using too large a ratio of <r d to cr in the cross- diffusion term adversely 
affects compressible bound ary- layer predictions. Aside from inclusion of the 
cross- diffusion term, the analysis is the same as that presented by Wilcox. 35 * 

I he log layer is the region sufficiently close to the solid boundary for neglect 
of convective terms and far enough distant for molecular diffusion terms to be 
dropped. Fn the log layer, the equations of motion based on the k-<j model 
simplify to the following. 


(u) Mach HU 


Figure A. I: Comparison of computed ami measured 
pressibic flat plate boundary layers; aj = .5; - 
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(A .2) 
(A- 3) 
(A 1) 


(A. 5) 


pT=p v ,r ui (A. 6) 

I hf quantity n r is friction velocity defined as y/r w /p v where r w is surface 
Rhear stress and p tt , is density at the surface. Also, T„ is surface temperature, 

ifl surface heat flux and C r is specific lieat at constant pressure. Finally y is 
difltance from the surfarr. 

We now change independent variables from y to w. Consequently, derivatives 
transform according to 

d du d d 

< A7 > 

With this change of dependent variables, we replace Kquations (A 3)-(A 5) 
hy the following. 


l /r T 2 J 

(A8) 

^ TT ^ * 2 4 — ^ 

du 2 

(A. 9) 

^ d 2 u; ^ <r d dkduf f u fip 2 kw 

du 2 F du du f f> F p*uj ' !! 

(A. 10) 

Integrating Equation (A. 8) yields the temperature, and hence (he density, as 
a function of velocity and Mach number based on friction velocity, M r = u r /u,„. 


£ = "f - [s &)'♦*£) + *• (*)) <*•"» 


Next, we assume a solution of tin* form: 

pk = l>,„u? 


(A. 12) 
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where 1 ’ is a constant to he determined. Substituting Equations (A ll) and 
(A. 12) into Equation ( A.9), and noting that Al? = 2 l'A/ r \ leads to the following 
cubic equation for V. 


F l> + (7 - l)/V r <T*A/^r] r 2 = I (A. 13) 

As can easily be verified, when <£ 1 , the asymptotic solution for 1' is 


r - 




(7-1)/ 


2 0* J Mr 


finally, in terms of l', Kqiiation (A. 10) simplifies to 
d J u> aj dp 


°-n - -f + I" - i"‘\ -TF = 0 

du 1 p du 1 1 p w uf I 


(A. 14) 


(A. 15) 


(xiiiibiiiing Lquations (A ll) and (A 12) yields the density as a function of 
velocity and 1 \ 


/V 

P 


ilzJi 

2 




1 + ( 7 - l)Pr,aTM? 


(A. 16) 


Equation (A. 16) assumes a more compact form if we introduce the freestream 
velocity, A bit more algebra yields 


Pw _ 1 d- It ii — A*jr* 

~~ 1 + ( 7 - l)/V/ <rTA/? 

where 

v = u/U^ 

A*' = ^y^/V/ 

It = -Pr r q u ,U 00 /(C p T w r Ul ) 


(A 17) 


(A. 18) 


Using Equations (A. Id), (A. 17) ami (A. 18), and retaining terms up to O(Atj), 
Equation (A. 15) assumes the following form, 


Z± \—!L <M f K'td’^hh) 2 1 

<r [l t Uv- A*v* du J [lAtfe-AVd '*' -0 

where the constant A' w is defined by 


Kl — “ 


(7- i)/Vi(:b>-/V^)<T- 
2 (T 


J Mr 


(A. 19) 


(A. 20) 
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and k is Karnian's constant. Hocanso V^/u, » I, we can use tin- VVKH method 
to solve Equation (A. 19). Noting that w decreases as v/Vnu increases, the 
asymptotic solution for u> is 


C 


+ If® - A V] ‘ ’ " 104,0 ,/4 exp [- K„ u* /»i r ] (A. 21) 


where (* is a constant of integration and u* is defined by 


«• i - / a 

— — = — sin I — 


2 A 2 v — 


Wo 


V/7T 


n - /y \ 
+ 4/1V 


(A. 22) 


Combining Equations (A. 2), (A 12) and (A. 21), we can relate velocity and 
distance from the surface. 


j[ I + Up- A , v'} < 7 *‘ f -" , \x V lK w u*/« T ) <h> 


<11 

re r .. 


(A 2;i) 


VVe integrate by parts to generate the asymptotic expansion of the integral in 
Equation (A. 23) as f/oo/u r — ► oo. Hence, 


[1 + 


2 21 (* + 2<Td/ff)/4 


eXp[/C«*/Mr) 


A’toiC'y 

F»r 


(A. 21) 


Finally, we set the constant of integration (' ~ V u * Taking the natural 
log of Equation (A. 24), we conclude that 


u 

Mr 


I 

K . 


-fn 


(?) + * 


(A. 25) 


The quantity is tlie effective “constant” in the law of the wall defined hy 

fl + 2<J w /»)/4 


b w = n + - 2 - 6 . f -C-) 

\/W 


(A. 26) 


where B is a true constant. 

As discussed by Wilcox, 29 the difference between A’ u . and I he Harman con- 
stant is of no great consequence. However, the variation of /.f w with the density 
ratio has a large effect on both skin friction and the predicted law of the wall. 
As clearly demonstrated by Wilcox, 11 the magnitude of the power to which 
p/p w is raised determines the degree of distortion of the velocity profile. For 
example, the Standard Jfc-f model is equivalent to a k-u> model with = 2<r. 
This corresponds to an exponent or 5/1, while the unmodified k-co model has 
an exponent of 1/4. As shown in Figure A.I, there is virtually no distortion of 
the compressible law of the wall when = (). Hy contrast, for the (r d = 0.3 
and cr<i = 0.5 cases we have used a ~ 0.0 and <r - 0.5, respectively, so that the 
corresponding exponents are 1/2 and 3/4, respectively. The figure shows in- 
creasing distortion of the compressible law of t he wall as <r (i increases. Although 
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not shown, the k-t model departs even farther from the compressible law of the 
wall. 

As a final comment, if we had written the cross-diffusion term as 


Cross diffusion = 


a a 0(pk) du> 
v dxj dxj 


(A 27) 


the exponent in Equation (A. *2(3) would be reduced to 1/4, independent of the 
value of a**. This is true only for Lite const ant- pressure case since pk is ap- 
proximately constant in the log layer. The same statement cannot be made 
for boundary layers with nonzero pressure gradient, however. A shortage of 
Contract funds precluded testing this alternative. 
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